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The vortex structure is studied in light of MgB2 theoretically based on a two-band super- 
conducting model by means of Bogoliubov-de Gennes framework. The field dependence of the 
electronic specific heat coefficient 'y{H) is focused. The exponent a in 'y(H) oc H°' is shown to 
become smaller by adjusting the gap ratio of the two gaps on the major and minor bands. The 
observed extremely small value a ~ 0.23 could be reasonable in this two-band model with the 
gap ratio ~ 0.3. 
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Much attention has been focused on the recently dis- 

I covered MgB2 because of its relatively high supercon- 
ducting transition temperature Tc 39K and simple 
crystalline structure. General consensus obtained so 
far is that the electron-phonon interaction is mainly re- 
sponsible for the pairing mechanism in this system be- 

' cause the large isotope effect is observed.^' '^^ There is, 

\ however, little consensus as to the microscopic descrip- 
tion for the record high due to the electron-phonon 

\ interaction. 

Apart from the much debated pairing mechanism, it is 

\ rather urgent to determine the precise pairing function or 
gap function realized in MgB2 . There are several impor- 

\ tant, but conflicting experimental data concerning the 
superconducting energy gap A, ranging from the strong 

\ electron-phonon coupling ~ 5 to extremely 

weak coupling value 2A/k-QTc ^ 2. These come from the 

\ earlier experiments, such as position-dependent tunnel- 
ing or Raman experiments.^^ More recent experiments 

\ show unequivocally that these two gap values come from 

• a single sample and converge to definite values: the 
\ larger and the smaller A5 whose ratio Ag/A/, falls 

around 0.3'^0.4. These experiments include photoemmi- 
; sion (AL=5.6meV, As=1.7meV, As/Al=0.30),^) the 

• T-dependent specific heat analysis (A5/Al=0.27)^' ''^ 
I and tunneling experiment (As/Al=0.42).^^ Through 

these analyses, they are able to obtain systematic and 
smooth T evolutions of each gap value. This implies that 
the two gap structure is an intrinsic property in MgB2. 

According to the band structure calculations,^' 
there are two distinctive Fermi surface sheets; one is a 
two-dimensional cylindrical Fermi surface arising from 
cr-orbitals due to Px and Py electrons of B atoms and 
the other is a Fermi surface coming from 7r-orbitals due 
to Pz electrons of B atoms. They are weakly hybridized 
with electron orbitals of Mg atoms. Since the cr-orbital is 
strongly coupled to the in-plane B-atom vibration with 



symmetry simply because the hopping integral be- 
tween the tJ-orbitals is modulated by this bond stretching 
motion. On the other hand, it is shown by the band cal- 
culation-'^''^ that the 7r-orbital is weakly coupled with this 
phonon mode. Thus it is quite conceivable that these two 
Fermi surfaces with different electronic characters have 
different energy-gap values if this particular in-plane vi- 
brational mode is responsible for the attractive interac- 
tion which induces superconductivity in MgB2. 

Here we are going to analyze the field dependence 
of the T-linear electronic specific-heat coefficient 7(H) 
in the superconducting mixed state by investigating the 
vortex lattice structure in two-band superconductors, ft 
is known that ^{H) is a sensitive and useful quantity to 
reflect the gap structure through the zero-energy exci- 
tation spectrum inside and outside the vortex core.^^"^^^ 
In particular, the exponent a in "f{H) oc iJ" at low fields 
reflects the nodal structure of the superconducting gap 
at the Fermi surface, playing a vital role to identify the 
gap function. Several recent specific heat experiments 
on MgB2 show a very small exponent^' '^^' '^^^ a ~ 0.23, 
implying that on increasing H the zero-energy density 
of states (DOS) in the mixed state quickly recovers its 
normal state value, compared with those in d-wave su- 
perconductors^ '^"'^^■' with a ^ 0.5 or clean limit s-wave 
case'^^' ^^^-^ with a ~ 0.7. Since there are no definitive 
reports which claim a line or point node in the supercon- 
ducting gap in MgB2, this small exponent remains mys- 
tery and requires a proper explanation by microscopic 
calculations. This is one of our purposes in this pa- 
per based on the microscopic theory of Bogoliubov-de 
Gennes (BdG) framework. '^''' 

We start with a model pairing Hamiltonian for a two- 
band superconductor described by tight binding form: 

H = Hq + Hpair, (1) 

Ho ^ {-ijjj - fJ"tSij)al^^ajaj, (2) 
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with the nearest neighbor (NN) hopping integral 

Uj.y = tjexp[i— A{r)-dr], (4) 

where A{r) is the vector potential and (f)o = hc/2e is 
the unit flux. The two-dimensional square lattice whose 
lattice constant is unity is assumed. The index 7 denotes 
the two bands j = L and S. Assuming the singlet pair- 
ing, we can derive the BdG equations for 7 = i and S 
in a standard way: 

(5) 

where 

The gap equation is given by 

= '^9jydY{ri) (7) 

7' 

with the order parameter 

= -X]^7E(''i)"7e(^i)tanh^. (8) 

e 

The local density of states (LDOS) at site i for the 7 
band is calculated by 

€ 

+ \v^,{r,)\''5iE + E^,)}. (9) 

We assume an isotropic s-wave pairing for both bands 
7 = L and S characterized by the order parameters (the 
energy gaps) (ii(Ai) and dsi^s)- The attractive in- 
teractions are chosen as gLL ^ 0, gLS = 9sl 7^ and 
9ss = 0, namely in eq. (5) the gap A5 on the S-band is 
induced by the Cooper pair timncling via 9ls- As for the 
normal state band parameters we take II =ts = t{= 1) 
and /XL = — 1 and us = +1, thus the Fermi surface for 
7 = L{S) is close (open) around the F-point. The DOS 
for both bands is same at the Fermi level. As two vortices 
are accommodated in a unit cell of Ng, x atomic sites, 
the applied magnetic field is given by Hjsf.^xNa = 2<i'o/^f ■ 
By introducing the quasi-momentum of the magnetic 
Bloch state we obtain the wave function under the peri- 
odic boundary condition for a large number of unit cells 
(See detailed numerical calculations in ref. 19). 

First, we study the case qll = 2.0 and gLS = 0.6, 
which gives Al = 0.322 and As = 0.086 at zero field. It 
is designed to adjust the gap ratio for MgB2(As/AL = 
0.27). If we consider a single band superconductor with 
the small gap A = 0.086, the superconductivity vanishes 
at the following magnetic field discussed below, since 



(a)-ff30x30 

\dL{r)\/\dL\H=o \ds(r)\/\ds\H^o 




Fig. 1. Spatial profiles of the order parameter |(ii(r)| for the L- 
band and \ds{r)\ for the 5- band at lower field -ff3ox3o(a) and 
at higher field -ff ig x 18 (b) • They are normalized by the zero field 
values |dz,|if=o(= 0.145) and \ds\H=o{= 0.055), respectively. 
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NL(r,Er^ 0)/Ar(EF) Ns(r, E ~ 0)/N(Ef) 
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Fig. 2. Zero energy local density of states Ari:,(r, E ~ 0) for the L- 
band and Ng (r, B ~ 0) for the 5-band at lower field H30 x 30 (a) 
and at higher field -ffisx 18(b). They are normalized by N{Ef), 
the normal state DOS at the Fermi level. 



H > Hc2- But, in this two band superconductor, the 
small gap superconductivity survives in the S-band be- 
cause of the cooper pair transfer g^s- 

The spatial profiles of the order parameters dL{r) and 
ds{r) arc shown in Fig. 1 where the unit cell of the 
square vortex lattice is displayed. Vortices are accom- 
modated at the center and four corners. It is seen that 
the vortex core radius for the L-band (S'-band) is small 
(large) and the depression of \d{r) \ is apparent along the 
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Fig. 3. Field dependence of 7(_ff) for As/Al = 0.27. Points of 
7(i?) (circles), 73 (H) (triangles) and 7^^ (if ) (squares) are numer- 
ical data. The thick line is fitting for lower field data of f{H). 
The thin line is fitting by 7(if) ~ 7Ar(i?/-ffc2)". In the low field 
the thick(thin) dotted line is fitting for ■ys{H){'yi^{H)). 



NN direction, in particular, for the 5'-band. By increas- 
ing H, \d{r)\ is further suppressed as is seen in Fig. 1 
where the core radius is widen. The suppression by H is 
eminent in the S'-band. 

The corresponding spatial profiles of the LDOS are 
shown in Fig. 2, where N^ir, E Q) and Ns{r, E Q) 
have a peak at the vortex center and the ridges con- 
necting the vortex cores are clearly seen. While the 
high density of states is concentrated at the vortex core 
in NL{r, E ^ 0), it rather spreads out in Ns{r, E ^ 0). 
This is because the vortex bound states are highly con- 
fined in the i-band vortex corresponding to the narrow 
core radius while in the S'-band vortex the core states are 
loosely bounded. The spatial profiles for NL{r,E ~ 0) 
and Ns{r,E ~ 0) are resemble to those of the low- 
field case and the high-field case in the single band su- 
perconductor (see Fig. 1 and Fig. 2 in ref. 13). In 
Ns{r, E ~ 0), the low energy states extending from vor- 
tex cores overlap with each other, and the LDOS is sup- 
pressed along the line connecting the NN or next NN 
vortices. With increasing H, the effect by the overlap 
becomes eminent, and the LDOS is reduced to the flat 
profile Ns{r, E 0) /N{Ef) ~ 0.5 in the S-band {N{Ef) 
is the total DOS in the normal state at the Fermi level). 

The spatial average of A^L(r-, - 0) and Ns{r,E '-^ 0) 
gives rise to the total DOS under a given field, leading 
to 7(i?) which is defined by 



with 



j{H)=jl{H)+js{H) (10) 



7L,S{H) = {NL,s{r, E ~ 0))r:unitcell. (H) 



We have done extensive computations for various 

HN.^xN-,,ciises. In Fig. 3, it is seen that 'y(H) is described 
by a power law: 7(i?) cx with small a. If only the 
low field points are fitted, we obtain a = 0.38(thick line). 
The fitting by j{H) ^ Jn{H/ Hc2)°' under the condition 
that 7(i?) is reduced to the normal state value gives 
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Fig. 4. Field dependence of "/(H) for Ag/Az, = 0.27(circles), 
0.49(triangles) and 0.997(squares). (a)Fitting lines for the low 
field data. (b)Fitting lines by 7(ff) ~ ^n{H/Hc)°'. Points 
of numerical data are the same in both figures. In the insets, 
As /Ai -dependence of a is shown for each fitting case. 



a = 0.33(thin line). The small exponents a, or the sharp 
rise of "fiH) in small fields, can be attributed to the S'- 
band contribution js{H) cx H°-^", while -/l{H) oc 
in the L-band. That is, the small a is due to the overlap 
of the low energy states outside of vortex cores at the 
S-band. Physically it is because the energy gap for the 
S-band is suppressed by a weak field, while the total su- 
perconductivity is maintained by the larger energy gap 
up to Hc2- This intuitively appealing picture is actually 
confirmed by the present microscopic calculation. This 
is, however, diflFerent from the two independent gaps with 
different transition temperatures and different Hc2- In 
such a case, we would have double transitions and j{H) 
would be a simple addition of two independent curves, 
which has a kink structure at the lower Hc2- This is not 
the case for MgB2. 

To study the dependence of a on the gap ratio As/Al, 
we perform the calculation for various pairing parame- 
ter sets. In Fig. 4, we show the As/Ai-dependence of 
7(iJ) behavior. There, we show the results using the two 
kinds of fitting; the fitting for low field data in Fig. 4(a), 
and overall fitting by -/{H) - Jn{H/Hc2)" in Fig. 4(b), 
while the numerical data are the same in both figures. 
In the insets, we show the Ag/Ai-dependence of a in 
each fitting case. In the hmiting case Ag/A^ ^ 1, a 
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is reduced to the exponent in the single band case. In 
Fig. 4(b), it gives a ^ 0.7 in accord with the previous 
quasi-classical calculation.^^) In both cases, a monoton- 
ically decreases with decreasing As/Al- It should be 
emphasized that wc may identify the gap ratio Ag/Ai^ 
uniquely by measuring the electronic specific heat under 
varying external field, providing a rather unique spectro- 
scopic method for determining the gap ratio. In Fig. 4, 
for the gap ratio As/Al ^ 0.3 observed by the several 
groups with different methods, wc obtain small exponent 
a, as in the specific heat data on MgB2. 

There are several factors which might alter our con- 
clusion on the relation a vs. As/A^. 

(1) Wo assimic that the DOS in the normal state for 
each Fermi siirfacc sheet is equal. According to the band 
calculation by Bclashchcnko et al?'^^ the ratio of the two 
DOS is 0.55 (TT-band) : 0.45 (cr-band). This small difl?cr- 
ence causes potentially to alter our conclusion, but not 
in an essential way. 

(2) It is assumed that in the minor S-band there is 
no direct attractive interaction gss = 0. The gap in the 
S-band is exclusively induced by the Cooper pair tun- 
neling process via g^s 0- According to Kortus et al,^^ 
the electron-phonon coupling in the S'-band due to the 
bond stretching mode is smaller than that in the L-band, 
but yet non-vanishing. Thus gss might not completely 
vanish in MgB2. We might regard it vanish as a first 
approximation because our conclusion relies exclusively 
on the gap ratio As/Al- The effect of gss will be 
studied in the future for more details. 

(3) We comment on the small discrepancy of the expo- 
nent a between our calculation based on BdG theory and 
that of the experiment. Since our parameters belong to 
a rather quantum limit case (the coherent length ^ is an 
order of the atomic lattice constant), the quasi-classical 
calculation^^' is more appropriate for MgB2. We be- 
lieve, however, that the overall relation a vs. As/A^ 
is not greatly altered in that calculation. We will study 
this case in a future publication. 

Haas and Maki^^^ analyze an anisotropic s-wave pair- 
ing state in connection with MgB2. Their single band 
model is designed to describe the anisotropic supercon- 
ducting properties such as the upper critical field or the 
penetration depth. According to the recent penetration 
depth measurement^^) for single crystals of MgB2, the 
anisotropy of the penetration depths ior H || c and H ± c 
in the hexagonal crystal is almost absent, which is at odd 
with the prediction by Haas and Maki. Since their sin- 
gle band model is similar to our two band model in the 
sense that the gap anisotropy is implemented in the re- 
ciprocal space in Haas and Maki or implemented in the 
energy space in ours. In order to fully describe the three 
dimensional superconducting nature in MgB2 our model 
should consider the anisotropic s-wave pairing state for 
both major and minor bands, which may better explain 
the above penetration depth experiment. 

We speculate that the present multi-gap model may 
have potentially wide applicability. It is quite usual that 
a superconductor has a multiple gap because the under- 
lying Fermi surface consists of multiple sheets, on each 
of which the gap value could be different. It is true even 



for elemental metals. MgB2 may belong to an extreme 
case. To reveal this feature, the measurement otj(H) is 
demonstrated to be a useful tool. The analysis for the 
p-wave pairing case,^^) focusing on Sr2Ru04, is reported 
in ref. 24. 

In conclusion, wc have evaluated the exponent a in the 
T linear specific heat coefficient jiH) cx H" for a simple 
two-band superconductor and succeeded in reproducing 
the extremely small a = 0.3, as in observed in MgB2, 
by taking the two gap ratio As/Al = 0.3, each coming 
from the different Fermi sheets. Thus we conclude that 
the gap functions are distinctly different for the different 
Fermi sheets, the major is the cr-band (px and py charac- 
ters of B atoms) while the minor is the 7r-band (pz char- 
acters of B atoms), yet each gap being isotropic on its 
own Fermi sheet. Thus we do not need exotic anisotropic 
gap function for describing superconductivity here. This 
two-band feature is intrinsic in MgB2. 
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